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1 Introduction 

We consider the motion of a passive scalar advected by a random velocity field V(t, x) = 
(Vi(t, x), • • • , Vdit, x)). The governing equation is 

^T = V(t,»(t)) (i) 

where V(£, cc) is a mean-zero, time-stationary, space-homogeneous random incompressible 
velocity field. 



In certain situations, it is believed that the convergence of the Taylor-Kubo formula (J14 
) given by 



oo 



J {E[V t (t, 0)^(0, 0)] + E[V t (t, 0)^(0, 0)]} dt (2) 

o 
is a criterion for convergence of passive scalar motion to Brownian motion in the long time 
limit. Indeed, it has been shown that the solution of 

^=I V (Jp **(«)), x e (0) = (3) 

converges in law, as e — > 0, to the Brownian motion with diffusion coefficients given by the 
Taylor-Kubo formula when the velocity field is sufficiently mixing in time (see [[F], ||, ||, 
0). Moreover, the solution of (|3|) converges to the same Brownian motion for a family of 
non-mixing Gaussian, Markovian flows with power-law spectra as long as the Taylor-Kubo 
formula converges (see 0]). In this paper, for the same family of power-law spectra, we show 
that, when the Taylor-Kubo formula diverges, the solution of the following equation 

^- = e^ s V(±x £ (t)), - £ (0) = 0, (4) 

with some 8 =£ 1 depending on the velocity spectrum, converges, as e — > 0, to a fractional 
Brownian motion (FBM), as introduced in |H| (see also ||13|| ). 



We define the family of velocity fields with power-law spectra as follows. Let (fl, V, P) 
be a probability space of which each element is a velocity field V(£, a;), (t,x) e R x R d 
satisfying the following properties. 

H 1) V(t,x) is time stationary, space-homogeneous and centered, i.e., E{V} = 0, and 
Gaussian. Here E stands for the expectation with respect to the probablity measure 
P. 

H 2) The two-point correlation tensor R = [R%j] is given by 

Rijit, x) = E [Vi(t, *)^(0, 0)] = J cos (k • a;)e" |k|2 ' 3 *R ij (k)cik (5) 

J R d 

with the spatial spectral density 

R ( k ) = | k |2a+d-2 ^ " ~^f) ' ( 6 ) 

where a : [0, +oo) — > R + is a compactly supported, continuous, nonnegative function. 
The factor I — k ® k/|k| 2 in (g[) is a result of incompressibility. 




<x = 



H 3) a < 1, (5 > and a + /3 > 1. 

The function exp (— |k| 2/3 t) in (^) is called the time correlation function of the flow V. 
For (5 > 0, the velocity field lacks the spectral gap and, thus, is not mixing in time. As the 
time correlation function is exponential, the Gaussian velocity field is Markovian in time. 

Because the function a has a compact support we may assume, without loss of generality, 
that V is jointly continuous in both (t,x) and is C°° in x almost surely. For a < 1, the 
spectral density R(k) is integrable in k and, thus, ©-© defines a random velocity field 
with a finite second moment. The exponent a is directly related to the decay exponent of 
R. Namely |R|(0, x) ~ |cc| a_1 for |cc| 3> 1. As a increases to one, the decay exponent of R 
decreases to zero. 

Our main result is summarized in the following theorem. 

Theorem 1 Under the assumptions H 1)- H 3), the solution of eq. (fjj with the scaling 
exponent 

5 .__ P 

a + 2/3 -I 

converges in law, as e tends to zero, to a fractional Brownian motion B//(t) that is to a 
Gaussian process with stationary increments whoe covariance is given by 



E [B H (t) ® B H (t)] = Dt : 



2H 



(7) 



with the coefficients D 



D 



Rd 



Hkp" _ 1 + | k |2/3 
1^120+4/3-1 



k®k\ a(0) 



|kl d - 



dk 



and the Hurst exponent H 



l/2<H=l/2+ a+ ^ - <1. (9) 



Remark. Molecular diffusion can be added to the equation of motion so that instead of (JTJ) 
we may consider an Ito stochastic differential equation 

dx(t) = V(t, x(t))dt + V2~KdB(t) 

with B(t), t > the standard Brownian motion, independent of V and k > 0. This however 
would not influence our results. 

2 Multiple stochastic integrals 

By the Spectral Theorem (see, e.g., []J) we assume without loss of any generality that there 
exist two independent, identically distributed, real vector valued, Gaussian spectral measures 
Vi(t, •), I = 0, 1 such that 

V(t,x)= fv Q (t,x,dk), (10) 



where 

V (t, x, dk) := c (k • x)V (t, dk) +c 1 {k- sc)Vi(t, dk) 

with co(0) = cos (0), ci(0) = sin (0). Define also 

Vi(£, x, dk) : = -c x (k • x)V (t, dk) + c (k • a?)Vi(t, dk). 

We have the relations 

dVo(t,x,dk)/dxj = kjVi(t,x,dk), (11) 

dVx^x.d^/dxj = -kjV {t,x,dk). (12) 

Clearly / Vi(t, x, dk) is a random field distributed identically to and independently of V. 
We define the multiple stochastic integral 

I ■ ■ ■ J i)(k u - ■ ■ ,k N )V h (t 1 ,x 1 ,dk l ) ®...<g> V lfr (t N ,x N ,dk N ) (13) 

for any Ii,---,In £ {0, 1} and a suitable family of functions i[) by using the Fubini theorem 
(see (pl|) below). For tpx,- ■ ■ , ip^ e S(R d ), the Schwartz space, and h,- ■ ■ ,In £ {0, 1} we set 

/ ••• / V ; i(ki)---^v(kAf)V il (ti,a3i,dki) g> • • • <g> Vi N (t N , x N , dk N ) (14) 

:= / V'lOfrOVjifajiCijdki) <g> •••® / V'Jv(kiv)V; JV (i J v,a; J v,dk i v). 

We then extend the definition of multiple integration to the closure 7i of the Schwartz space 
S((R d ) N ,R) under the norm 

:=|---|v(ki,---,MV(k' 1 ,---,k^ v ) (15) 



E 



Vj^ti, Xi, dki) ® • • • ® Vi N (t N , x N , dk N ) ■ V^(ti, x x , dki) ® • ■ • ® Vi N (t N , a?7v, dk' 



iV/ 



The expectation is to be calculated by the formal rule 



E [V hi (t, x, dk)V Vtil (t', x', dk')] = e" |k|2/3|t -*' l 5 z/ c (k ■ (x - x'))R itl ,(k)5(k - k')dkdk'. 

This approach to spectral integration follows 

When i = (i 1; • • ■ , id), h, ■ ■ ■ , id G {1, 2, ■ • • , d} is fixed and 1 = (Zi, ■ • • , Zat), Zi, • • - , Zjv £ 
{0, 1} we shall denote the corresponding component of the stochastic integral by ^y. 

Note that ^y G H N (V) - the Hilbert space obtained as a completion of the space of 
iV-th degree polynomials in variables / ip(k) V(£, x, k) with respect to the standard L? norm. 

Proposition 1 For any (t 1 ,x 1 ),- ■ ■ ,(t N ,x N ) G R x R d and p > 0, ^i belongs to L p (Vt) 
and 

(E|^ M n 1 / p <C(E|* 1 , i | 2 ) 1/2 (16) 

with the constant C depending only on p, N and the dimension d. Moreover, ^i is differ- 
entiable in the mean square sense with 

W M (£i, ■■-,t N ,x 1 ,---,x N ) = (-1) 1 ' J ■ ■ ■ J k^ki, • • • , k N ) (17) 

V hjil (tx, xi, rfki) • • • Vi^^.(^-, Xj, dkj) • • • Vi N}iN (t N , x N , dk N ). 

The proof of Proposition ^ is standard and follows directly from the well known hyper- 
contractivity property for Gaussian measures (see, e.g., ||, Theorem 5.1. and its corollaries), 
so we do not repeat it here. 

The field V is Markovian i.e. 



E 



JiP(k)V l (t,x,dk)\V- 00 , s ~\ =|e-l k l 2 ^ t - s V(k)VKs,a ; ,rfk), 1 = 0,1, 



for all ip G S(R d , R), where V a ,b denotes the a-algebra generated by random variables V(t, cc), 
for t G [a, b] and x G R d . 

To calculate a mathematical expectation of multiple product of Gaussian random vari- 
ables, it is convenient to use a graphical representation, borrowed from quantum field theory. 
We refer to, e.g., Glimm and Jaffe || and Janson H. A Feynman diagram T (of order n > 
and rank r > 0) is a graph consisting of a set -B(^ r ) of n vertices and a set E(F) of r edges 
without common endpoints. So there are r pairs of vertices, each joined by an edge, and 
n — 2r unpaired vertices, called free vertices. B(F) is a set of positive integers. An edge 
whose endpoints are m, n G B is represented by mh (unless otherwise specified, we always 
assume m < n); and an edge includes its endpoints. A diagram T is said to be based on 
B[T\ Denote the set of free vertices by A[T\ so A[T) = T \ E^). The diagram is com- 
plete if A(F) is empty and incomplete, otherwise. Denote by G{B) the set of all diagrams 
based on B, by Q C {B) the set of all complete diagrams based on B and by Qi(B) the set 
of all incomplete diagrams based on B. A diagram T' G Q C (B) is called a completion of 
T G Qi{B) if E{T) C E{T'). 



Let B = {1, 2, 3, ..., n}. Denote by T\k the sub-diagram of J-, based on {1, • • • , k}. Define 
Ajt(jF) = y4(jF| fc ). A special class of diagrams, denoted by Q S (B), plays an important role in 
the subsequent analysis: a diagram T of order n belongs to Q S (B) if A fc (jF) is not empty for 
all k = 1, ...,n. 

We shall adopt the following multiindex notation. For any P £ Z + , multiindex n = 
(jii, ■ ■ ■ ,rip), |n| stands for J2 n p- If -P' < -P we denote nip/ := (ni, • • • ,np>). In addition if A; 
is any number we set n • k := (m, • • • , np, A;). 

We work out the conditional expectation for multiple spectral integrals using the Markov 
property (|T8|). 



Proposition 2 For any function ip (zTi and h,- ■ ■ ,In £ {0> 1}? *i; ' ' ' > ^n £ {!>■''> ^}> 



E 



^(ki, • • • , k N )Vi uh (t, xi, dki) • • • t^ W)ijv (t, x N , dk N ) | V_ 



(19) 



^ ■■■ expl- Yl \k m \ 2(3 {t- s)\ifj(k 1 ,---,k N )V StXl ,-,x N {dk 1 ,---,dk N ;J : ') 

F&g({l,...,N}) J J [ m&A{T) J 

V S)Xl ,-,x N (dk ir -- } dk N ;J r ):= ]J Vi mAm (s } x m} dk m ) (20) 



x n 

mn£E(T) 



1-e 



-([kml^+lknl^)^-*) 



E [^ m ,i m (s, #m, dk m ) V^ n (s, a; n , dk„)]. 



Proof. Without loss of generality we consider ^(ki, 



^N, 



l j4l (ki)---l Ajv (k JV ) for 



some Borel sets Ax, ■ ■ • , A N . 

JNTote that \i{t,Ai) = V°(t,A) + V\{t,Ai) where Vf(t, ■) is the orthogonal projection 
of Vi(t, ■) on L 2 _ oot and V[(t, ■) its complement. Here L 2 6 denotes L 2 closure of the linear 
span over V(s, cc), a < s < 6, a; e P d . The conditional expectation in ( |T9"[ ) equals 



T&g{{l,...,N}) f^h£E{F) m£A(F) 

The statement follows upon the application of the relations 

V°(t,A) = /V |k|2/3( '- s) V,(s,dk) 

A 

and 



E 



V^^OV^S) 



J J 6i, v {E [Vi(t, dk) <g> V r (t, dk')] - E [v°(£, dk) ® V°(£, dk'j\ } 



A B 



3 Proof of tightness. 

We begin with the following lemma which shows, among other things, that the family of 
continuous trajectory processes x e (t), t > is tight. 

Lemma 1 For the family of trajectories given by (fQ) we have 



limE \{x E {t) - x e {r)) ® (x £ (t) - x £ (r))} = D|£ - r 

e|0 



2H 



where H , D are given by ^), (||) respectively. 

Proof. Thanks to the stationarity of the path x £ (t) it is enough to prove the lemma for 
r = 0. By the stationarity of V(s,ex(s)) (fTl|]), we write 

t 
\im~E[x £ (t)(g)x £ (t)] = \ime 2 [ ds /"e [V(s', ex{s')) ® V(0, 0)] ds' (21) 



which equals 

JV 

2^2i n + n N (22) 



n=l 

where 

t 

~^Z s s„_ 1 

l n = e n+1 J ds J ds x --- J E[W„_i(5i,---,s ft ,0)®V(0,0)]ds n 
00 

and 

W (si,aj) = V(si,aj) 

W n (si, • • ■ , s n+ x, x) = V(s n +i, x) ■ VW„_i(si, ■■■ ,s n ,x) for n = 1, 2, • • • 

with the remainder term 

t 

K N = 2e N+2 J dsj d Sl --- / E [W N ( Sl , • ■ ■ , sjv+i, eaj(s^+i)) ® V(0, 0)] ds N+1 . (23) 
00 

Estimates ofX n . Elementary calculations show that 

limJi = Dt 2H , for a + (3 > 1. (24) 

Since V is Gaussian we have that 



We now show that 



EX ra = 0, for even n. 



limEJ n = 0, for odd n. (25) 

ej.0 



Set 

E s „ +1 Wn-i.^Si, • • • , s n , x) := E W n -i,i(si, • • • , s„, cc) I V- 



oo,s n+1 



The i, j-th entry of the matrix X n is given by 

t 

~^5S s Sn-1 

e n+1 J dsj d Sl --- J E[E W n _ ltl { Sl ,- ■■ ,3^0)^(0, 0)]ds n . (26) 



The conditional expectation in (|26| ) can be expressed in terms of spectral measures of the 
velocity field. To do so we introduce first the so-called proper functions of order n, a : 
{1, • • • ,n} — > {0, 1} that appear in the statement of the next lemma. The proper function 
of order 1 is unique and is given by cr(l) = 0. Any proper function, a', of order n + 1 is 
generated from a proper function a of order n as follows. For some p < n, 



o'(n + l) 
a'{k) 
o-Xp) 





cr(k) for k < n and k ^ p 

l-a(p). (27) 



In other words, each proper function a of order n generates n different proper functions of 
order n + 1. Thus, the total number of proper functions of order n is (n — 1)!. In the sequel, 
we sometimes write a^ instead of a(k). 

Lemma 2 Let n > 1 and Si > s 2 > ■ ■ ■ > s n > s n+ i, i 6 {1, • • • , d}, x e R d . We have then 
that 

E, B+1 W n _i ii (si, ■•-,s n ,x)= (28) 

H / ■" /^i>( k i 5 "--,kn)exp{- Y, |k m | 2/3 (sn-s n+ i)}x 
P n -i{F)Q{F) II Vi m ,<r m (s n +i 1 x,dk m ), 

meA n {F) 

where ip\ n J are some functions with sup \<p\ n J\ < 1 

Pn-i(F) = li[ E |k m |)exp{- E \k m \^( Sj - Sj+1 )} (29) 

3=1 \m€Aj(F) j mZAjiT) 

Q(f) = _Ji E p m ,. m (o,rfk m )v v;iCTm; (o,rfk m o] • (30) 

mra'e£(f) 

The summation is over all multiindices i of length n, whose first component equals i, all 
T ' E Q s and all proper functions a of order n. 



Before proving Lemma |2|, we apply it to show (|25|) . By ( p8|) 

t 

~^25 s s n -i 

I ds J dsi • • • / E W„_i,<(si, • • • , s n , 0)ds„ 







t 



for i = 1, • • • , d. Here, adopting the convention s n+ i := 0, we set 

^(ki,---,k 



5 ■ lv ray 



fdsf- J (pf a (ki,-- -,k n ) x n exp<^ - E |k m | 2/3 (sj - s i+ i) > dsi • • • ds n 

3=1 [ meA,-(:F) J 



s n 



/ dsi • • • / n exp < - E |k m | 2/ S' >ds x --- ds n 

j=l [ meAj(T) ) 

and 

„_if/ \ l-exp{- E |k m | 2 M 

^-^) = n e mU -^fL_ — 

It is elementary to check that, due to \(PiJ\ < 1, 

I^SCki,--',^)! < 1. (31) 

By Lemma [| the left hand side of (^) equals 



71* 



#J(k 



2e^ E / ^ / • • • / > V k ' W K-i(F)Q(F)V 



meA„(7) 



II V im , am (0,dK 

mgA n (Jf)U{n+l} 



(32) 
Here the summation extends over all multiindices i = (i\, ■ ■ ■ , i n +x) such that i\ = i, i n+ \ = j, 
all Feynman diagrams T <^Q S and all proper functions a of order n. Note that 

1 - e~t/ £2S C 

— r- £ ^ < 33 > 

for all positive e, £. Here and in the sequel C stands for a generic constant independent of 
e. C in fl3~3|) is also independent of £ > 0. Thus, the absolute value of (|32]) is bounded by 



. n+l-2<5\T^ /" /" Pn-l,e(^J tt o(fe m k m i)dk m dk m i 

^ J "J p25, V ^ 1J - fc 20 - 1 l j 

^ £ £ + 2^ K m ' — -, ^m 

with 

meA,^) 



Pn-l^C? 7 ) : = n 

Using the fact that 



U 25 + E fern"' 



E "'m A , 

<G^i%, Vm^MT), (35) 



P-25 _i_ v k 2/3 ~ F 2S 4- k 2f3 
m&Aj{T) 



9 



we have 



Pn-lA?) <C\\ £ ~" +k 



e 2S + E km 

m£A n (T) 






Vm^e^iJ). 



(36) 



Bounds (|36|) and (0) imply that 



\? n \ < Cte 



n+1-26 



e n 



S \1n 

* [eV + h, 



TeQc 



LLrC'r 



mm'£E(F) 



(> + k 2Jiy™+8 m , mn k 2a-l 



(37) 



where q m are certain nonnegative exponents satisfying 

E 9m = n - 1 

mm'eE^) 

and & m ,m n = if m ^ m„ and = 1 otherwise. 

The integrals appearing in the expression ( |37"D are of the form 



(38) 



, s \1 

* [ee+k 



dk 



( £ 25 + k 2py+r k 2 



q+r U2a-1 



(39) 



for some q > and r G {0, 1}. They may diverge or remain bounded as e J, depending on 
the exponents q, r. If q, r are such that the integral diverges then 2f3(q + r) + 2a > 2 + q 
and, consequently, 

(1 + k) 11 dk 

x — — - < +oo. 



;i + my+ r k 2a ~ i 

u 

In either case, the integral ( |39| ) is bounded from above by (7£ n ( 2 -"-2/3)/(o!+2/3-i) f or ^ < « + 2/? 
after a change of variable k'- = kj/e^ in case (^) diverges as e — > 0. Therefore 



_n+l-2<5 



n(2-q-2)3) 



-1 



XJ < Cte n+i -"e -+2/5-1 _i < Cte^w^r 



(40) 



which vanish as e j for n > 3. 

Estimates oflZjq. By 



727 s 



■SAT 



7^ = 2e iV+2 J dsj d Sl ---JE [E SN+1 W N ( Sl , ■ ■ • , Sjv+1 , escCsiv+i)) ® V(0, 0) 



By the Cauchy-Schwartz inequality we get that 

[^ivl 2 < 4t 2 e 4 C 1 -«) +2JV 'E|V(0,0)| 2 x 



ds 



N+l- 



(41) 



max E 

0<s<t/e 2S 



■■ • E SJV+1 Wiv(si, ■ ■■ ,s N , sn+i, ex(s N+1 ))dsi ■ ■ ■ ds N+1 



!>S1>--->SJV+1>0 



10 



The stationarity of the Lagrangian velocity field implies that the maximum in ( fH"|) is equal 
to 



max E 

0<s<t/e 2S 



8 

/ ds' ■ ■ ■ E Wat(si, • ■ ■ , s^, 0, 0)dsi ■ ■ ■ dsN 



s'> Sl >->s N >0 



< 



(42) 



C max E 

0<s<t/e 2S 



8 

ds' /•••/ B VW N - 1 (si,-",s N -. 1 ,SN,O)dsi'--ds N 



xE|V(0,0)| 2 . 

s'>si>--->sjv>0 

Here the hypercontractive property of the Gaussian measure is used. Subsequent applications 
of Lemma |2] to ( f42|) yields the upper bound 

2 
t 2 



C ^s E 



E 



</v(k l5 ■ • • , k N )P N (F)Q(F) J] ^wJO, dk, 

m£A N (T) 



(43) 



with some IV'i.o-I < 1- The summation above extends over all Feynman diagrams T G Q s , the 
relevant proper functions a and multiindices i. 
Thus, we have 



-^'-^ ram' 



o\k rn k ra i)dk rn (ik ri 



(44) 



Here the summation extends over all possible completions of with JF e Q S ({1,- ■ ■ ,N}), 
T' G G S ({N + 1, ■ • ■ , 2iV}). The product is over all edges of any completion of T U T' . 
Arguing as for ( p7[ ) we obtain that 



\Tl N \ 2 < Ct 4 e 2N+ ^- 2 ^ Ell/ 
for some q m > with 



£0 + k m \ dkr, 



■p <j:t — — ~-- " 
~* ^ mm' 

E q™ = 2N - 



r 25 i h?P I U2a-1 



(45) 



(46) 



Moreover, 



\n N \ 2 < Ct'e^^-^e- 2 - 1 ^^ 1 < Ct'e^h-^ l - 25 \ 



(47) 

which vanishes as e j. for a sufficiently large N. In conclusion, we proved that the left hand 
side of (|2l|) tends to Dt 2H as e j 0, provided that a + /3 > 1 (see (||)). 

By the hypercontractivity property of the L p norms over Gaussian measures we also know 
that for any p > 1 and T > there exists a constant C > 



E\x e (t)-x £ (s)\ p <C(t-s) 2Hp 



(48) 



for any T > t > s > 0, e > 0. 

Proof of Lemma pi. We prove the lemma by induction. The case n = 1 is obvious 
by choosing ip\ = 1. Suppose that the result holds for n. For the sake of convenience we 
assume with no loss of any generality that s n+ 2 = 0, then 



E W n+hi (s ll • • ■ , s„ + i, x) = E |V(s„ + i, x) ■ VE Sn+1 W n>i (s 1 , ■ ■ -, s n , x)j 



(49) 



11 



By virtue of the inductive assumption we can represent E Sn+1 W U)i using (]28f) and as a result 
(|49|) becomes 



EEo 



J ■■■ J ^ n j(k l5 • • • , k n ) exp{- J2 |k m | 2/3 (s„ - s n+1 )}P n ^(F)Q{F) (50) 



V (s n+ i, x, dk n+1 ) ■ V \ ]J ^ m ,«T m ( s »+ii ^, dk TO ) 



To calculate (|50D we decompose each V^ ;i (s, cc, <ik) as 

V ff)i (s, x, dk) = V°i(s, x, dk) + V^(s, x, dk) 

where 

V^(s, x, dk) = e-M 20 ^V a ,(t, x, dk) 

is the orthogonal projection of V a ^ on V-oo,t- Expression (|50D becomes 



(51) 
(52) 



EEc 



J ■■■ J <p$(k u ■ • • , k n ) exp{- £ |k m | 2/3 (^„ - Sn+OJ^-i^g^/C^) 



with 



/C(.F) := 



mei„(f) 



m£A„(J) 



(53) 



£'={ej} 
jeA n (^)U{n+l} 



The term corresponding to Qj = 1 vanishes, as is clear from the following calculation, 

E{J ■■■ J <p < $(k l ,---,k n )Pn-i(F)Q(F) (54) 

Vl(s, x, dk n+1 ) ■ V J] Kw™0> *> rfk « 

\meA„(.F) 



VE 



/ • • • / ^(k 1; • • • , k n )Pn-i(^)g(-^)Vi( S , as, dk n+1 ) I] %.*>, ^ 5 d ^) 1 

meA„(f) J 



by homogeneity of the velocity field. By (Q)-(Q) 



V (s,x,dk n+1 ) ■ V< J] ^^(SjXjdkn. 

E k m' • V (s, x, dkn +1 ) Y[ V <T rn> jim (s,x,dk r , 



(55) 



m'eA„(f) 



where 



m£i„(f) 



1 — aw if m! = m 



a. 



o,, 



otherwise. 



(56) 
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By (|52|) , (|5l|), (|55|) and the definition fl2"9|) , (|53|) further reduces to 



E 



(n) ™m' ,i 






(57) 



exp{- ]T \k m \ 2(3 s n+1 }P n (F)Q(F) [] K^ ti Jt,x,dk n 

m&A(F) meA{F) 



n [i- e - 

pqeE{T>) 



dkpi^+mi 2 ")^-*) 



E 



\ m/ , v (0,0,rfk p )\ m „ 9 (0,0,rfk, 



with d\^ m > = and c^+i,™/ = aj 1 and all incomplete Feynman diagrams T' based on the set 
A n {T) U {n + 1}. 

Lemma 2 follows with 



¥>{;+>!, • • • , k n+1 ) : = ^g(k 1; ...,k„) " m/ ^;, n I 1 






m'£A n (F) 



pqeE{T>) 



4 Proof of weak convergence 

It is easy to see that the Gaussian processes 



y e {t) :=e f V(s,0)ds t>0. 



(58) 



converge weakly to the fractional Brownian Motion B#(£), t > given by (|7|). In addition 
we have 

limsupE|y e (t)| p < +oo 

for any p > 1, t > 0. 

We now prove that 

limE {[x Ml (ti) - x £jil (t 2 )] Pl • • ■ K, M (£ M ) - x e , iM (t M+ i)] PJVI } 

e|0 



E {[%(*!) - SH,i 1 (t 2 )] P1 ■ • • [^mM - S ff>iM (t M+1 )]^} 



(59) 



which, in conjuction with the tightness, identifies the fractional Brownian motion Buit) as 
the limit. Equation ([59]) is a consequence of 

lim |E{[x ejil (t!) - x £>h (t 2 )] pi ■ ■ ■ [x £ , iM (t M ) - x £ , iM (t M +i)] PM - (60) 

sj.0 

[y s ,M - y s ,M\ pl ■ ■ ■ [ye,i M it M ) - y e , iM (t M+ i)] PM }\ = o 

with y £ {t) = (y e ,i(t), ■ ■ ■ ,y £t d(t)), which, in turn, follows from the next lemma. 
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Lemma 3 For any positive integers M , pj, multiindices ij G {1, • • • , d} Pj with j = 1, • • • , M 
we have 



lim 

ej.0 



E 



'(pi) i 



t(pm). 



r( P l) 



(pm)i 



Zfflfa h)- Z2£(0, t M ) - W^(t 2 , h) ■ ■ ■ W™(0, t M ) = 0. (61) 



Here for any integer N > 1, multiindex i = (ii, • • • , in) £ {!•••> ^} ant ^ t > s we define 

r r N 
Z { E f(s,t) := e N J ■ ■ ■ J J] ^.(sp^x^p))^! • • -dsjv 



Ajv(s,t) 



and 



W. 



r r N 
g°(M) :=e N J ■■■J H V lp (s p ,0)d Sl ■ --ds N , 



Ajy(s,t) 
with A N (s,t) := {(si, • • • , s N ) : t/e 2S > s x > ■ ■ ■ > s N > s/e 25 }. 

Proof. To avoid cumbersome expressions that may obscure the essence of the proof we 
consider only the special case of M = 1 and t\ = t, t 2 = 0. The general case follows from 
exactly the same argument. We shall proceed with the induction argument on p\ = P. The 
case when P = 1 is trivial because the stationarity of the relevant processes implies that the 
expression under the limit in fl61|) vanishes. By (|4q) we know that 

\imsupE\Z^)(0,t)\ q < +oo, Vg > 1. 

ej.0 



Suppose that ( Ell ) is true for P > 2 and that 



limsupE|^~ 1) (0,*)| ? <+oo, Vg > 1. 

ej.0 



Like (|2~2T) we write 



with 



N-l 



EZ { £ y(0,t)= J2Zn(0,t)+K N (0,t) 



(62) 



(63) 



n=0 



J n (0,t) := e p+n+1 / • • -/E | E S2 W^(si n) , £: r( S2 )) f[ K p (s p ,ex(s p )) I ds£ n) ds 2 • • -rf Sp (64) 
A(™)m^ y p=2 > 



A^'(0,t) 



K N (0, t) := e p+N+1 / • • 7 E I E <i,w< ( s ?°» ^ wO) ft ^( s p> ^( s p)) \ ds[ N) ds 2 ■ ■ ■ ds P . 



A£"(0,t) 



(65) 



Here 



A?(M) := {(s ( r\s2,---,s P ) : t/e 25 > s£ n > > s 2 • ■ • > Sp > s/e 25 } 

with Si := (si,i, • • ■ , Si, n +i)- We write £ > Si > s, if t > Si > s n > s, where s = (si, • • • , s r 
is any ordered n— tuple in the sense that s± > ■ ■ ■ > s n . 
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The argument of the proof of Lemma |2| and (62) imply that 



and 



limJ n (0,t) = 0, 7i >1 

ej.0 



\imTZ N (0,t) = 

e|0 



for N sufficiently large. Thus EZ^; (0, t) has the same limit as the term 

J (0,t) :=e p+1 I '•■■ fE\v h (s 1 ,ex(s 2 ))f[V ip (s p ,ex(s p ))\ds 1 ---ds P . 



(66) 



A<?>((M) 



For (|5T)| ) we use a generalization of the argument of the proof of Lemma [2| Let us introduce 
some additional notation. For any multiindex i = (ii, • • • , i p ) and p > 1 we define W[' n by 
induction as follows. We set 

W£°... tip (si, ■■■,s p ,x) := Vi^sx, x) ■ ■ ■ V ip (s p , x) - E{V h (s!, x) ■ ■ ■ V ip (s p , x)} 

and 

W^+l( Sl , ■ • • , Sp _ 1? s p n+1 \ x) := VW^.. )ip ( Sl , •••,s("U). V(s p , n+2 , x) 

for any ordered (n + l)-tuple s p n > = (s P) x, • • • , s p>n+ i) < s p _i and (n + 2)-tuple s p n+1 *> = 
(s Pt i, • • • , s Pin+ i, s P)n+2 ). Expanding the left hand side of (|6|) like (|^) we obtain that 

J (0,t) = 
e p+1 /■•• fE{y il (s 1 ,CT(s 2 ))K i2 (s 2 ,ea;(s 2 ))}E [] V ip (s p , sx{s p )) \ d Sl ds 2 - ■ -ds P + 

\(o),A U= 3 J 



A™(0,t) 



jV-1 



^X lin (0,t)+^ li7 v(0,t) 

n=0 

where 

Xi, n (0,t):= (67) 

e p +" +1 / • • •/ E | E 83 W?%( Sl , s { 2 n \ex(s 3 )) n V^(s p ,eaj(s p )) 1 d Sl ds { 2 n) ds 2 ■■■ds P 

^i,jv(0,t):= (68) 

1 J ■■■ J E \ E s 2tN+1 WiJ 2 (si,s {N \ex(s 2tN+1 )) Y[V ip {s p , ex(s p )) \ ds t dsi N) ds 3 - ■ -ds P , 

A^ N \0,t) { P=3 J 

Ap n) (0,t) := {( Sl , S { 2 n \s 3 , ■■■,s P ): t/e 2S > Sl > s 2 n) > ■ ■ ■ > s P > 0}. 

We represent the conditional expectations appearing in (|6?D and fl6"8"|) using a generalization 
(Lemma |j) of Lemma 0. 
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£ p+iV - 



To formulate it we need a generalized notion of a proper function, which we call a p- 
proper function. Let p be a positive integer. The p-proper function of order 1 is unique and 
is given by a(i) = 0, i = 1, • • ■ ,p. Any p-proper function, a', of order n + 1 is generated from 
a p-proper function a of order n as follows. For some q < p + n, 



cr'(p + n + 1) 
a'{k) 

At) 



o 

cr(k) for k < n + p and k ^ q 

l-a(q). (69) 



We also distinguish a special class of Feynman diagrams Qg{B) : a diagram JF of order 
n + p belongs to G^(B) if A^{T) is not empty for all k — p, ..., n +p. 

Lemma 4 For any positive integer p, s\ > • • - > s p _i > sj^ 1 ) > s, a multiindex i = 
(ii, • • • , i p ) G {!,••-, d} p we /iat»e 



E.sW( 



p,n—l 



[s 1 r--,s p ^- 1 \x) = j:J---J<p^ ) (k 1 ,...,k p+n ) (70) 



ex p{- z) ik m i 2/3 (s p , n - s)}p Pi „_i(^ r )Q(^ r ) n ^.^(s.^^y. 

mgyt„ + p(f) m£A„ +p (f) 

where ^pf„ are functions satisfying \^pf„\ < 1 and 

n+p-l / \ _ 

^p,n(y J = _[_[ I 2^ l^ml I ex Pl~ 2^ |^m| [ s p,j-p ~ s p,j-p+l)fi 
j=p \meA,-(.F) J m&AjiF) 

Q(?) = _II E [^ m , CTm (o,rfk m )^ m/i(Tm ,(o,rfk m o] . (n) 

mm'eB(f) 

TTie summation is over all multiindices j = (ji, ■ • • ,jn+p)> such that ju, = i, a// T $zQ v s and 
all p-proper functions a of order n. Here by a convention s Pt o := Sp_i. 

The proof of Lemma f| is exactly the same as that of Lemma ^| and is omitted. 

Continue the proof of Lemma 3 using Lemma 4 we have that X (0, t) is asymptotically 
equal to EZ^j (0, £) and 

e P+l J ••■ J E< K il (si,£:^(s3))V r j 2 (s2,£:a[;(s3)) JJ V^Sp.ea^Sp)) > ds 1 ds 2 • ■ • ds P 
A P (0,i) I p= 3 J 

is asymptotically equal to EZJ; (0, £), as e j 0. Repeating the above argument p-times 
we obtain (|6T|). Finally the hypercontractivity properties of the L p norms over Gaussian 
measure space imply that (^) holds with P — 1 replaced by P n 
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